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Abstract — The purpose of this paper is to explain the phenomenon of symmetry breaking for optimal 
functions in functional inequalities by the numerical computations of some well chosen solutions of 
the corresponding Euler-Lagrange equations. For many of those inequalities it was believed that the 
only source of symmetry breaking would be the instability of the symmetric optimizer in the class of 
all admissible functions. But recently, it was shown by an indirect argument that for some Caffarelli- 
Kohn-Nirenberg inequalities this conjecture was not true. In order to understand this new symmetry 
breaking mechanism we have computed the branch of minimal solutions for a simple problem. A 
reparametrization of this branch allows us to build a scenario for the new phenomenon of symmetry 
breaking. The computations have been performed using Freefem++. 
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1. Symmetry breaking for Caffarelli-Kohn-Nirenberg inequalities 

In this paper we are interested in understanding the symmetry breaking phenomenon 
for the extremals of a family of interpolation inequalities that have been established 
by Caffarelli, Kohn and Nirenberg in [1|. More precisely, let d ^ 1, # € (0, 1] and 
define 

r^, a iP — 2 d — 2 , d , 2d 

&{p,d) :=d— — , a c := — — , b = a-a c + -, and p (&,d) := - — — . 

2p 2 p d — 2v 
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Notice that 2* = p*(l,d) = if d ^ 3, while we set 2* = °° if d = 1, 2. For any 
<i ^ 3, we have 

(K©(p,d)<0<l 2^p^p*{p,d) s£2*. 

We shall assume that # £ (0, 1], a < a c , and 2 ^ p ^ p*(p,d) if either <i ^ 3, or 
J = 2 and # < 1, or J = 1 and ■& < 1/2. Otherwise, we simply assume 2 ^ < oo; 
& 6 [a,a + 1] if <i ^ 3, b G (a,a + 1] if d = 2, and b G (a + ^,a + \] if d = I . Under 
these assumptions, the Caffarelli-Kohn-Nirenberg inequality amounts to 



AT 



bp J IS^-iKp- 2 )/? V JM rf kl 2 " J \J^ d M 2(c,+1 



(1.1) 

with A = (a — a 2 ), for all functions w in the space obtained by completion of the set 
^(W 1 \ {0}) of smooth functions with compact support contained in U. d \ {0}, with 
respect to the norm 



II |xr a Vw||? 2 ^ + ||bcH a+1 



L 2( R d) -r || |*| ' 'W|| L 2 

We assume that Kckn(# 5 P>A) is the best constant in the above inequality. We will 
denote by K^ KN (#,p,A) the best constant when the inequality is restricted to the 
set of radially symmetric functions. 

According to [2], the Caffarelli-Kohn-Nirenberg inequality on W 1 can be rewrit- 
ten in cylindrical variables using the Emden-Fowler transformation 

s = log|xi, co = j-r G S^ 1 , u(s,(o) = \x\ ac ~ a w(x). 

and is then equivalent to the following Gagliardo-Nirenberg-Sobolev inequality on 
the cylinder ^:=Rx S rf_1 , namely 

II m IIlp(^) ^ K C kn(^,,P,A) [||VM||^ (y) +A||M|^2 (y) J ||m||l2(^ (I- 2 ) 

for any u G H x (^€\ Here we adopt the convention that the measure on S rf_1 is the 
uniform probability measure. 

The parameters a <a c and A > are in one-to-one correspondence and we have 
chosen to make the constants Kckn and K£ KN depend on A rather than on a because 
in the sequel of the paper we will work on the cylinder 'if. Also, instead of working 
with the parameters a and b, in the sequel we will work with the parameters A and p. 

Note that u is an extremal for (11.2b if and only if it is a minimizer for the energy 
functional 
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Radial symmetry of optimal functions in dl.lb . or symmetry to make it short, 
means that there are optimal functions in (11.21 which only depend on s. Equivalently, 
this means that Kckn(#,A,/?) = K£ KN (#, A,p). On the opposite, we shall say that 
there is symmetry breaking if and only if Kckn(#,A,/?) > Kq KN (-&,A,p). Notice 
that on the cylinder the symmetric case of the inequality is equivalent to the one- 
dimensional Gagliardo-Nirenberg-Sobolev inequality 

INIlp(r) ^ k ckn(#> a >.p) ( II v "IIl 2 (m) + a II"IIl 2 (r)) IHIl 2 m 



for any u £ H (K). The optimal constant K GKN (#, A,/?) is explicit: see Q. 

Symmetry breaking of course makes sense only if d ^ 2 and we will assume it 
is the case from now on. Let us summarize known results. Let 

, nN d-l (2#-l)p + 2 , , . (N-l)(6-p) 
toW.-* I and A, W:= L_ii_W. 



Symmetry breaking occurs for any A > Aps according to H4I31 (also see |2J for pre- 
vious results and JH if d = 2 and ■& = 1). This symmetry breaking is a straightfor- 
ward consequence of the fact that for A > Aps, the symmetric extremals are saddle 
points of the energy functional Q l K in the whole space, and thus cannot be even local 
minima. 

If ■& = 1, from @ , we know that symmetry holds for any A ^ A+(p) . Moreover, 
according to $7), there is a continuous curve p h-> A s (p) with lim p ^2 + A s (p) = °° and 
A s (p) > a\ for any p G (2,2*) such that symmetry holds for any A ^ A s and there is 
symmetry breaking if A > A s . Additionally, we have that lim p ^2* A s (p) ^ a c if d ^ 
3 and, if d = 2, lim^oo A s (p) = and lim p ^ 00 /7 2 A s (/j) = 4. The existence of this 
function A s has been proven by an indirect way, and it is not explicitly known. It has 
been a long-standing question to decide whether the curves p — > A s (/?) and the curve 
p — > Aps (p, 1 ) coincide or not. This is still an open question. Let us notice that for all 
p e (2,2*), A*(p) ^ A s (p) ^ Aps(p, 1) and that the difference A F s(p, 1) — A*(p) is 
small, and actually smaller and smaller, when the dimension d increases. For more 
details see O. 

According to 0, existence of an optimal function is granted for any # € 
(Q(p,d),l), but for # = Q(p,d), only if K C kn(#,A,/>) > K GN (/?), where K GN (/?) 
is the optimal constant in the Gagliardo-Nirenberg-Sobolev inequality 



2 K GN (fl,p,A) 2 &(p,d) „ ,,2(i-e(p,rf)) v uiw, 

\\ U \\lJ>(RN) ^ |g^-l|(p-2)/p H V "Hl 2 (R") II M Hl 2 (R'V) VMetl (K j. 

A sufficient condition can be deduced, by comparison with radially symmetric func- 
tions, namely K^ KN (&(p,d), A,p) > Kqn{p), which can be rephrased in terms of A 
as A < A GN (p) for some non-explicit (but easy to compute numerically) func- 
tion p i — y A GN (p). When # = &(p,d) and A > A GN (p), extremal functions (if 
they exist) cannot be radially symmetric and in the asymptotic regime p — > 2 + , 
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this condition is weaker than A > K-ps{p,®{p,d)). One can indeed prove that 
lim p _ > .2 + Afs(P)©(p,^)) > hm / ,^2 + Aq N (/?). Actually, for every # in an interval 
(0, i9#), &d < 1, and for p £ (2,2 + e) for some e > 0, sufficiently small, one can 
prove that Aps(p,&(p,d)) > Ag N (p); see (9l for more detailed statements. Hence, 
for # £ (®(p,d), 1), close enough to &(p,d) and p — 2>0, small, optimal functions 
exist and are not radially symmetric if A > Aq N (/>), which is again a less restrict- 
ive condition than A > Aps(p, #)• See @ f° r proofs and ifTOl for a more complete 
overview of known results. 

The above results show that for some values of # and p there is symmetry break- 
ing outside the zone, parametrized by A, where the radial extremals are unstable. 
So, in those zones the asymmetric extremals are apart from the radial extremals and 
symmetiy breaking does not appeal - as a consequence of an instability mechanism. 
What is then the mechanism which makes the extremals lose their symmetry ? The 
goal of this paper is to understand what is going on. A plausible scenario is provided 
by the numerical computations that we present in this paper. They have been done 
using Freefem++. 

The paper is organized as follows. In Section [2] we expose the theoretical setup 
of our numerical computations. In Section [3] we describe the algorithm. Section |4]is 
devoted to the numerical results and their consequences. 

Our numerical method takes full advantage of the theoretical setup and provides 
a scenario which accounts for all known results, including the existence of non- 
symmetric extremal functions in ranges of the parameters for which symmetric crit- 
ical points are locally stable. Although we cannot be sure that computed solutions 
are global optimal functions, we are able to present a convincing explanation of how 
symmetry breaking occurs. 



2. Theoretical setup and reparametrization of the problem for # < 1 
with the problem corresponding to # = 1 

Let us start this section with the case # = 1 and consider the solutions to 

-Au + pu = u p ~ 1 in . (2.1) 

Any solution u of (12.11 ) is a critical point of Q 1 ^ with critical value Q^u] = 
Up to multiplication by a constant, it is also a solution to 

-Au + ptu ^-u p ~ 1 =0 in < T, (2.2) 



where 
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Since Q 1 ^ [A u] = Q 1 ^ [u] for any A 6 R \ {0}, (12.11 ) and (12.21) are equivalent. For sim- 
plicity, we will consider primarily the solutions of (12. It . 

Let us denote by h^* the unique positive symmetric solution of (12.11 ) which 
achieves its maximum at s = 0. We know from previous papers (see for instance [0 
and references therein) that the positive solution is uniquely defined up to transla- 
tions. As a consequence it is a minimizer of Q 1 ^ among symmetric functions and 

G^iKu,*] = ll M Ai,*llLP(f(f) = 1/Kckn(1' A 4 '/')- 

Our first goal is to study the bifurcation of non-symmetric solutions of (12.11 ) 
from the branch (w^,*)^ of symmetric ones. Let f\ be an eigenf unction of the 
Laplace-Beltrami operator on the sphere S d ~ l corresponding to the eigenvalue d — 1 
and consider the Schrodinger operator ffl := — -jj^ + fj, + d — 1 — (j>— 1) u p ^ 2 whose 
lowest eigenvalue is given by X\(jl) = d- l+n-\np 2 (see ifTTH. p. 74, for more 
details). As in [4], let jUps be such that X\ (jUfs) = 0, that is 

d-l 



Mfs =4 



P 2 



We look for a local minimizer of by expanding = w^* + £(pf\ +o(e) 
in terms of £, for /I in a neighborhood of (jUfs)+- We find that QMm/x,* + £f\ <p] ~ 
£ 2 ((p,J4? (P)l 2 (V) as e — > 0. The problem will studied with more details in lTT2l . 

It is widely believed that for # = 1 the extremals for the Caffarelli-Kohn- 
Nirenberg inequalities (11.21 ) are either the symmetric solutions w^* for }l ^ /ips 
or the solutions belonging to the branch which bifurcates from u^s '■= W/x F s,*- Nu- 
merically, we are going to see that this is a convincing scenario. 

For # < 1, the Euler-Lagrange equation for the critical points of can be 
written as 



;i-0)*[v]+a 



#||v" p - 2 



in tf, (2.3) 



L"(%f) 



l|Vv||^ (y) ||Vv||^ (y) +A||v||^ 
where f [vj := M 2 — - — and x = jp-^ — = <2a[ 



MIl 2 ^) H V llL"( f #) 

Any solution u of (12.21 ) is also a solution of (12.31 ) if A = # — (1 — ■&) t [u]. Symmet- 
ric solutions w^* give rise to a symmetric branch of solutions jj. i-)- v A €>^ # = w^* 



of solutions for (12.3b . where 

A?(m):=0M-(1-0Mk/i,*]- 
A branch /j. \— > of solutions for (12.11 ) indexed by pL, normalized by the condition 
\\ufi WiJpg) = x can be seen as a branch of solutions of (12.21 ) and also provides a 



branch jU H> Vw#(u) = of solutions for (12.31 ) with 

A*0i):=d|i-(1-*)*["M]- 
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If we can prove that pi^-u^ bifurcates from \i \-t m^,* at jU = /ips, then we also have 
found a branch A h-> va of solutions of (12.31) which bifurcates from A h-> va,* at 



Ap* s := (hfs ) = * Mfs - (1 - #) < [«fs] 

as has already been noticed in Q- 

So, from the branch of solutions to <\2.2\ which bifurcates from wps we will 
construct a branch of solutions to (I2.3t that contains candidates to be extremals for 
the Caffarelli-Kohn-Nirenberg inequalities for # < 1 . Of course, nothing guarantees 
that the extremals for the inequalities (11.21 ) with # < 1 lie in this branch but, as we 
shall see, such a scenario accounts for all theoretical results which have been ob- 
tained up to now. In this paper we numerically compute the first branch bifurcating 
from the symmetric branch at jUps , then transform it to a branch for (12.31 and study 
the value of /^(ju) := Q^Jw/i] in terms of A^/i). 

An important ingredient in our algorithm is related to the fact that minimizing 
the first eigenvalues of Schrodinger operators —A — V under some integral constraint 
on V is equivalent to solving (12.ll ). One can indeed prove that 



_ inf inf f v \Vu\ 2 dy-f v V\u\ 2 dy 

l|v|| ww =i kgWAM h\u\ 2 dy 

with q = has a minimizer, and that the operator —A — V has a lowest negative ei- 
genvalue, — /i, with associated eigenfunction «„. Up to multiplication by a constant, 
V = u P p L /||w | l^Zgy so that solves (12.11 ); see for details, and also M13I14II . 



3. The algorithm 

Based on the theoretical setup of the previous section, we are now ready to introduce 
the algorithm used for the computation of the branch of non-symmetric solutions 
of (11.21 ) which bifurcates from /ips for # = 1 . 

1) Initialization: obtaining one point on the non-symmetric branch. The function 
u^^ is a saddle point of Q 1 ^ for any pt > /ips- Choose then /lo larger than jUps, but 
not too large, e > small and w a direction of descent for Q]^ at u^. Starting from 
Wjuo,* +£vv, we use the conjugate gradient algorithm to decrease the energy and 
search for a quasi-local minimum of Ql . The limit u^ solves (12.21 ) for some xo = 

<2L["/io]- Since its energy is lower than <2L[ m /jo,*]> it is certainly a non-symmetric 
critical point. 

2) A fixed point method. Critical points of Q 1 ^ can be characterized as fixed points 
of an ad hoc map as follows. 

1. Choose k > 0, p € (2,2*), q = and start with some potential Vb normal- 
ized by the condition: ||Vo||l*(#) = 1- 
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2. For any i ^ 1, define 

A,-(x) := inf ( / \Vu\ 2 dy — x / V/_i |m| 2 ] , 

II m IIl 2 (^) = i 

and get a minimizer U[ € H 1 ^) such that H^Hl 2 ^) = 1- 

3. Define V,-:= | 1 ^ 2 / 1 1 1 1 l^^t ) and iterate, by computing A, + i(x) from Step 2. 

The sequence is monotone non-increasing. Indeed, for any i ^ 1, we have 

A;+i(x) = inf ( / \Vu\ 2 dy — x / V;|w| 2 dy) 

II^IIl 2 ^) = i 

$5 / |Vw,| 2 cfy — x \ Vi\ui\ 2 dy 

inf ([ \V Ui \ 2 dy-x f V \uj\ 2 dy) 
V G L<?(<«f) V-/* •/«" / 

II^IIl*^) = i 



^ / \Vui\ 2 dy — x / V,_i \ui\ 2 dy = A,-(x) 

TV H 



The sequence is bounded from below, as an easy consequence of Holder's 

inequality: 

/ \Vu\ 2 dy — >c Vi\u\ 2 dy^ / |Vm| 2 dy — x ||V||l?(<^) , 

J'io 

and the r.h.s. itself is bounded by the Caffareli-Kohn-Nirenberg inequality: 

inf ( / |Vm| 2 dy — x\\u\\ 2 Plcn I =: ju(x) 
«eH'(?) \ J * / 

II m IIl 2 (<^) = i 

if we assume that ||V||l»(^) = 1. Indeed this amounts to 

II^IIl^+mMIMIl 2 ^) > *\\ u \\lp(<e) VmgH 1 ^), 

which is exactly equivalent to (11.11 ) up to a reparametrization of /I in terms of x. 
This scheme is converging towards a solution of 

— Am + ju(x)m = xVu , V = x\\u\\ 2 jj^,u p ~ 2 in"?f . 
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See O for details and H13I141 for earlier references. Notice that if we start with 
the potential Vq = u p ^ 2 /\\u^ found at the end of the initialization of our al- 
gorithm, then we find that /lo = jU(xo) and (as well as Vq) is a fixed point of our 
algorithm, such that x = Q^u^] < Q^u^,*]. 

The above iterative algorithm is a local version of a Roothan algorithm to com- 
pute a fixed point. We have run it using Freefem++ in a self-adaptive way, where at 
every step the computing mesh is based on the level lines of the previously computed 
function. This is important for large values of fi because solutions asymptotically 
tend to concentrate at some point. 

3) Building the branch in terms of x, starting from xq. We adopt a perturbative 
approach by modifying the value of the parameter x and reapplying the above fixed- 
point algorithm. 

In practice we take x = x$ — 77 for x\ > small, Vq = \u^ 1 /:> 2 /||«po HlpC*?)- ^ e 
get a new critical point with pt = p.{x). If f] has been chosen small enough, 
we still have that <2^[ M ,u] < Q]i[up.,*\- By iterating this method as long as pL > /i F s, 
we obtain a discretized branch of numerical solutions pt 1— > of (12.21) contain- 
ing UfjQ. Numerically, we check that Qjjw/t] < Gu[ M /i,*]> hence proving that is 
non-symmetric as long as /I > /ips, and such that converges to mfs as M tends 
to jUps- For simplicity, we adopt the following convention: we extend the branch to 
any value of /i > but observe that coincides with w M * for any /i < /Xps- 

We can do the same in the other direction and start with k = xq + 1] for 77 > 
small, take again Vq as initial potential, and then iterate, with no limitation on x. 
Again we check that QUw^] < for the discrete version of the branch cor- 

responding tO /I = jli(x), X > 

Altogether, we have approximated a branch which bifurcates from the symmet- 
ric one at /i = jUps- This branch is a very good candidate to be the branch of the 
global extremals for the Caffarelli-Kohn-Nirenberg inequalities for ■& = 1. The main 
reason for this belief is that if we start our algorithm for a jx close enough to /ips , 
we always hit the first (in terms of pt) branch bifurcating from the symmetric one, 
that is, the one bifurcating from /ips and observe that its energy is below the energy 
of corresponding symmetric solutions. On the other hand, the asymptotic value as 
H — > 00 is the one predicted by Catrina and Wang in |f2] for optimal functions. The 
estimates found in O indicate that the set of parameters in which a different branch 
of optimal functions would co-exist with the branch of critical points we have com- 
puted is remarkably narrow, and close in energy with the one we have found at least 
for /i = jUps- The existence of another, distinct branch of non-symmetric solutions 
which does not bifurcate from the symmetric ones, but still has the same asymptot- 
ics as ji, — s> 00, is therefore very unlikely. Hence we expect that our method provides a 
complete answer for optimal functions and for the value of the best constant in dl.ll > 
for = 1. 

Once we have constructed a discretized branch of solutions for # = 1, we use 
the transformation described in Section [2] to get a discretized branch of solutions 
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for (12.3I ). Same comments apply as for the case # = 1 (see lfl2l for the asymptotics 
as }i — > oo) and we expect that the computed solution is the optimal one for (ll.lt 
with ■& < 1. 

Let us finally notice that by lfT5l (also see lfl6l ). the extremals of (11.21 ) enjoy 
some minimal symmetry properties. They depend only on s and on one angle, which 
can be chosen as the azimuthal angle on S . For any d ^ 2, the problem we have 
to solve is actually two-dimensional, which greatly simplifies the computations. 

4. Numerical results 

Solutions have been computed for various values of p and # € [&(p,d), 1] in the 
typical case d = 5. We adopt the convention that solutions with lowest energy are 
represented by plain curves while symmetric ones, when they differ, are represented 
by dashed curves. Darker parts of the plots correspond to minimizers for fixed A, at 
least among the branches we have computed. Hence, define the functions 



where is the solution of (12.11) which was obtained by the method of Section[3]and 
Ufi^ is the symmetric solution of d2.ll ). For /I < jitps, these two functions coincide 
and their value is a good candidate for determining the best constant in (ll.ll ). At 
H = /Xps . non-symmetric solutions bifurcate from symmetric ones and for /I > /^ps 
the corresponding guess for optimal constants is given respectively by /* and jf. 

With the notations of Section 12 we know that A' ?=1 (/i) = }X. For # < 1, the 
parameter we are interested in is A and we look for the solution of (12.31 ) which 
minimizes Q^. The reparametrization of Section |2] and the bifurcation at /I = jltps 
suggest to still parametrize the set of solutions by /I and consider 



However, we will see that there is no reason why pL i— > A^(jU) should be monoton- 
ically increasing for jU > /ips ; this is indeed not the case for certain values of d, p 
and #. 

To illustrate these preliminary remarks, assume first that d = 5 and p = 2.8. For 
# = 1, the bifurcation of non-radial solutions from the symmetric ones occurs for 
H = n FS « 4.1667. For ■& ranging between 0(2.8,5) w 0.714286 and 1, branches 
can be computed using the reparametrization of Section |2] See Fig. 1. Although it 
is hard to see it on Fig. 1, right, the curves jj, h> (A 1? (/i),7' > ( i u)) may have a self 
intersection. We are now going to investigate this issue in greater details. 

In the critical case & = &(p,d), the limiting value J„ := lim fi -> 00 J & ( p ' tl > (n) 
along the branch of non-symmetric solutions corresponds to the best constant in 
Gagliardo-Nirenberg inequalities: 




[« M ] and jf{pL):=Q 




^(a^),/^)) and (A*>U*(M)) 



= k 



inf 

»eH 1 (M'')\{0} 



Jfjd |Vm| dx + J R d \u\ 2 dx 



2 



(J Rd \u\Pdx)i> 
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J*(/x)and7f(/x) 




Figure 1. Left- Plot of jX (plain curve) and jl i — >• ^ (jU ) (dashed curve) for d = 5, 

p = 2.8, # = 1. The branch of non- symmetric functions bifurcates from the branch of sym- 
metric ones for J l (jips) « 4.17 and 2^ FS [«fs] ~ 15.65. Right.- Plots of fi i— > J®(n) and 
H^-Jf(n) (dashed curve) ford = 5, p = 2.8, for # = 0.72, 0.75, 0.8, 0.85, 0.9, 0.95, 1. 



where k := (st}^ (1 - tf) 1 ^) \$ =& ( p j), according to 02). When $ = &(p,d), at 
least in cases we have computed, there are optimal functions (see Fig. 2.) if and 
only if 

A SC A GN := sup{Af {p4) (ji) : jf (p4) (ji) < /.} . 



7 # (,u) and/* (/j) 






/ I (f»Fs)=4(l"FS 






J " _— — / 






A A 9 


A 



: A 


A* 
' l FS 


















A'fs 





Figure 2. Case a" = 5, p = 2.8, # = 0(5,2.8). Left.- P/of o/ ju ^ (A°(ll),J*(n)) a/id 
(dashed curve) p: M- (A; (fl),jf (jtl))- Right-- Reparametrization p: Af (ju) (dashed 
curve) and p: h-> A tf (/i): f/iey differ for pi > pl FS « 4.17. 
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Fix some $o € (0, 1). Now we consider the subcritical regime, that is when p 
varies in the range (2,p*($o,d)). When p is close enough to p*($o,d), the branch 
stills bifurcates towards the left from the symmetric branch at Ap S = A d (jUps) and 
then, for larger values of Li, turns towards the right, crosses the symmetric branch 
(in the plane (A, [it]), not in the functional space), and then stays under this sym- 
metric branch for any larger value of Li. In other words, the map Li h-» A' 5o ( i u) is 
monotone decreasing in (jUfsjJUo) for some Llo > Mfs and then monotone increasing 
in (jUo,°°). See Fig. 3. This case is interesting, because when the non-symmetric 




2.845 2. 850 2.B55 2.860 2.B65 ~ I I — — — — < T — — 1 — — „ — — — — — i — — — — — f 

Figure 3. Plots for d = 5, p = 2.78, # = 0(5,2.8). Left- Plot of fi ^ (A*(n),J° (p.)). 



Right.- Reparametrization p h-> Af (ju) (dashed curve) and p i— > A^(p): the function Af 
is increasing for p < /ips, while the function A^ is decreasing on (pps,p\) and increasing 
forp>p h 

branch crosses the symmetric one, if the critical points are actually the optimal 
functions for dl.lb . then a symmetric extremal and a non-symmetric one coexist. 
Denote by Ai the corresponding value of A. In this case, for A < Aj the extremals 
of (11.2I) are all symmetric. At A = Ai there is coexistence of a symmetric and a 
non-symmetric extremal, and for any A > Aj the extremals are non-symmetric. Let 
us denote by the minimizer of Q^: the map A h-> is not continuous, since at 
A = Ai there is a jump. 

In terms of p, there is some /ij* G (0, jUfs) such that Af (juj" ) = Ai, and extremal 
functions for A < A\ are parametrized by p \-± (Af (ju),m^ *). There is also some 
jUi € (/A), 00 ) such that A f, (/ii) = Ai, and extremal functions for A > A\ are para- 
metrized by p (A' ? ( j u),m iU ). At A = A], it^* * and are two different optimal 
functions. 

In Figs. 4 and 5, we show the very different shapes of the symmetric and the 
non-symmetric extremals at A = Ai. 
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[O.ji/2] 



Figure 4. Case d = 5, p = 2.78, # = 0(5,2.8) and A\ = Af (/Ij*), corresponding to 
the crossing of the curve jX t— > (Af (jx),J^ (fx)) with the non-symmetric curve jX >— > 
(A^(jx),J^(fx)). Plot of the symmetric solution, u^**: Left - level lines, Right - 3d plot. 

JO,7t/2) 



\«1 




[O.ir/2] 



Figure 5. Case d = 5, p = 2.78, # = ©(5,2.8) and Ai = A^(fX\), corresponding to 
the crossing of the curve jX t— > (Af (jx),J^ (fx)) with the non-symmetric curve jX M> 
(A tf (n),J® (fx)). Plot of the symmetric solution, u^: Left.- level lines, Right.- 3d plot. 



When p takes smaller values in the range (2,p*(#o,c?)), the branch bifurcates 
towards the right and stays under the symmetric branch. In other words, the map 
jX i-> A^(ju) is monotone increasing for /I > /ips- See Fig. 6. 




Figure 6. Plots for d = 5, p = 2.7, # = 0(5,2.8). Left.- Plot of H ^ (A® (li),J 6 (jx)) 
and (dashed curve) jX i— > (Af (jtt),7*(ju)). Right.- Reparametrization jX ^Af(fx) (dashed 
curve) and jX h- > A^(jx). 



Symmetry breaking in Caffarelli-Kohn-Nirenberg inequalities 
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We may now come back to Fig. 1, right. An enlargement of the curves /I 1-4 
(A' ? (/i),/ !? ( i u)) clearly shows how one moves from the limiting pattern of Fig. 2 
(left) to the generic case of Fig. 3 (left) and finally to the regime of Fig. 6 (left) when 
■& varies in [®(jj,d), 1): see Fig. 7. 




Figure 7. Case d = 5, p = 2.8: curves ju i-> (A* (ju),/ 1 * (ju)) /or # = 0(2.8,5) « 0.7143 
(a), t> w 0.7213 (b) and t> « 0.7283 (c). 



Concluding remarks 

In this paper, we have observed that any critical point of is also a critical point 
of and can be rewritten as a solution of (12.11 ) up to a multiplication by a constant. 
Using reparametrizations, it is therefore obvious that extremals for (II. lb belong to a 
union of branches that can all be parametrized by /I. We have found no evidence for 
other branches than the ones made of symmetric solutions and of the non-symmetric 
ones that bifurcate from the symmetric solutions (when the number of eigenvalues 
of the operator M 1 changes as /I increases). Among non-symmetric branches, the 
first one is the best candidate for extremal functions in (11.1b . 

At this point we have no reason to discard the possibility of secondaiy bifurc- 
ations. Branches of solutions which do not bifurcate from the symmetric ones may 
also exist. Among the various branches, we have no theoretical reason to decide 
which one minimizes the energy, and the minimum may jump from one to another. 
This is indeed the phenomenon we have observed for instance in Fig. 3. 

However, the branch that we have computed is a good candidate for minimizing 
the energy. It is the natural one when one starts with small values of \x and tries to 
optimize locally the energy functional, and it has the correct behavior as pt — > °°. 
Known estimates, like the ones of O, show that there is not much space for unex- 
pected solutions in the range of parameters or of the energies. It is therefore quite 
reasonable to conjecture that the solutions that we have computed are the actual ex- 
tremals for Caffarelli-Kohn-Nirenberg inequalities and provide a complete scenario 
for the symmetry breaking phenomenon of the extremals, even if a complete proof 
is still missing. 
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